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Abstract 

We give an exposition of a technique, based on the Zwanzig pro- 
jection formalism, to construct the evolution equation for the reduced 
density matrix corresponding to the n-particle sector of a field theory. 
We consider the case of a scalar field with a gcj) 3 interaction as an ex- 
ample and construct the master equation at the lowest non-zero order 
in perturbation theory. 

1 Motivation 

The Hilbert space of a quantum field theory is the Fock space T constructed 
from the n-particle subspaces H n as T = 0^L O H n - ^ * s a wen known 
feature that the interaction terms cause mixing of the states between these 
subspaces. 

It is the aim of this letter to offer a treatment of the subspaces H n as 
corresponding to an open quantum system, that is we want to obtain an evo- 
lution equation for the reduced density matrix corresponding to the n-particle 
subspace, while treating the rest of the Fock space as an "environment". Our 
master equation is expected to contain (in general non-local in time) terms 
describing dissipation and noise. In addition, the fact that at any time the 
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system might leave the relevant subspace the density matrix will not be nor- 
malized to one. The quantity 1 — Trp re i(t) will correspond to the probability 
that at time t the n-particle subspace is "empty" . 

This is a problem of interest mainly for two reasons. First, in an attempt 
to explain the manifestation of particles as spatially localized excitations 
it is important to know the degree of validity of the particle description. 
Our approach of treating the n-particle sector as an open system provides 
measures of the modification in the equations of motion (dissipation) as well 
as of unpredictability (noise) inherent in such a coarse grained description. 

Our second motivation is the necessity of having a technique reduing 
the description of a field system at the particle level. This is important in 
the context of studying open field-theoretic quantum systems. We can derive 
master equations describing the evolution of reduced density matrix of a field 
system, by tracing out over the states of the environment. The environment 
can for instance be a thermal bath or another field. We are then naturally 
led to the question of what the evolution according to the master equation 
will imply for the particles corresponding to the field. Our preoccupation 
is particularly with the case of QED, where we can obtain the evolution 
equation for the spinor field by tracing out the states of the photon field [[]]. 
The translation of this into a particle description is related to the notion of 
environment induced decoherence [Q] and the possibility of dynamicalorigin 
of the charge superselection rules 0. 

In this letter, we are primarily interested in the mathematical aspects of 
the reduction of the field description into the particle level. We are therefore 
going to consider the simple case of a g<p 3 scalar field theory, even though 
we do not expect to obtain results of any physical interest. We should stress 
though, that the techniqua developed here is easily generalized to apply for 
any field theory that can be treated perturbatively and for all (even non- 
unitary) dynamics. 

2 Preliminaries 

The Zwanzig technique For the construction of the master equation we 
employ the Zwanzig projection technique ( for review and discussion see 
[0- |], |7|]). The main concept in this formalism is the representation of the 
coarse graining operation by an indempotent mapping P in the space of 
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states 

P^Prel = Pp (2.1) 

The irrelevant part of the state is then given by 

p irr = (1 - P)p (2.2) 
Considering unitary dynamics for the full theory 

t^ = Lp=[H,p] (2.3) 

we can show J3J that the evolution equation for the relevant part is given by 

i dp I edt)_ = pLprd{t) + PLe -(i-P)Lt pjrr(0) _ • f dTG ( r)prel{ t _ T ) (2.4) 
at Jo 

Here G stands for the kernel 

G(r) = PL(1-P)e iLr (l-P)LP (2.5) 

In equation (3.) the third term corresponds to non-local in time noise and 
dissipation while the second is essentially a time dependent "external force 
" acting on the relevant part of the state. 

In order to apply the Zwanzig method to our problem we need consider 
coarse-graining operations of the form 

Prel = P P = P nP P n (2.6) 

where P n is the projection operator corresponding to the n-particle subspace. 
The form of these operators is particularly transparent in the Bargmann 
representation, in which we are going to carry our calculations. We here 
make a small digression for purposes of establishing notation and assembling 
a number of useful formulas f| . 

The Bargmann representation In the Bargmann representation one 
can associate to each operator A on T a functional of a(x) a*(x) where 
a(x),a*(x) are elements of H (complex valued functions on a Cauchy hyper- 
surface Sigma of Minkowski spacetime). We actually have two choices for 
these functionals 

A{a\ a) = J2 ( n!r ^ni/2 /Ann(xi,...,x m I yi,...,y„) 

a*( Xl ) . . . a*(x m )a(y n ) . . . a( yi )d m xd n y (2.7) 
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in terms of the distributions A mn of the matrix representation of the operator, 
and 



A(a*, a) =J2 J Ann(xi, . . . , x m | yi, . . . , y n ) 
a*(xi) . . . a*(x m )a(y n ) . . . a(y 1 )d m xd n y 



(2.8) 



in terms of the distributions K mn appearing in the normal representation of 
the operator A 



A = J2 J K mn(*i, • • • , x n |yi, . . . , y m ) 
a f ( Xl ) . . . a t (x m )a(y 1 ) . . . a(y m )d m xd n y 

where a(x) and a^(x) satisfy the commutation relations 

[a(x),a t (x)] =5(x-x') 

The two functional are related by 



A(ct*, a) = A(a*, a) exp 



a* (yi)a(x)dx 



(2.9) 



(2.10) 



(2.11) 



The functional corresponding to the product of operators C = AB is given 
by 

C(a*,a) = J D(3*D(3A{a\(5)B{(3\a)e-^ p (2.12) 

In terms of the creation and annihilation operators a(x) and a^(x) the Hamil- 
tonian for a free field reads 



with 



Hq = i J dxdx'a\yi)h(x., x')a(x') 

/i(x,x') = y ^ e - ik(x ~ x,) cj k 



(2.13) 
(2.14) 



Here dk stands for the measure d 3 k/(27r) 3 . 

The evolution operator Uo(t) = e~ ltH ° reads in functional form 



Uo(a*, a; t) = exp 



J cfo;a*(x)A(x - x';t)a(x') 



(2.15) 



where 

A(x - x'; t) = J dke-^-^e-' 1 ^ 1 (2.16) 
We also note the fundamental forula of Gaussian functional integration 

DP* W(xi) . . . /3*(x m )/3(x' 1 ) . . . P(x' n )e-^ + ^ + ^ 

Xm+n 

-eT' (2.17) 



6f(x 1 )...5f(x m )6f(yf 1 )...5f(y^) 

Index notation It is more convenient to represent the functional using an 
index notation. To any function or distribution assign an abstract index to 
each of its arguments. The index is lower or upper according to whether the 
corresponding argument is integrated out with an a or an a* respectively. 
Hence we represent a(x) with a a and a*(x) with al and A(x — x'; t) with 
A b a . For example, the potential operator to encounter in the following 

V =: J dx^(j) 3 : (2.18) 

can be represented as 

V(a* , a) = (v a b c o: a a b a c + 3a a V a bc a b a c + 3a* a a* b V ab c a c + a a a b a* c V abc ^j 

(2.19) 

with the correspondence 

Vabc - I nLi ^ 77 ^e-( klxl + k2X2+k ^^(27r) 3 5(k 1 + k 2 + k 3 ) (2.20) 

V a bc ^/ni l7 ^e-(- k ^+ k ^+ k ^)(27r) 3 5(k 1 + k 2 + k3) (2.21) 

V ab c ~> / nLi (2 ^ ) i/ 2 e' t( ~ klXl ~ k2X2+k3X3) (27r) 3 ^(k 1 + k 2 + k 3 ) (2.22) 

Vabc ~> I nf=i ^ T7? e l ( klX1+k2X2+k « X3 )(27r) 3 5(k 1 + k 2 + k 3 ) (2.23) 

The inversion of the indices essentially amounts to the complex conjugation 
of the distribution. 
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3 The master equation 



Now in the Bargmann representation the projection operator P n corresponds 
to the functional 

P n {a\a) = -(a*a) n (3.1) 
n\ 

Substituting this into the pre-master equation(3.8) we obtain 



dt 

rt 

-I 



[H ,p(t)] + [P n VP n ,p(t)] 

fdr \-A n (r)p(t - r)Bt(r) - B n (r)p(t - r)Al(r) 
Jo 1 



+C n (r)p(t - r)C n (-r) + &Jr)p{t - r)Ct(-r)} + F res {t) (3.2) 

where the operator-valued kernels A, B>, C read 

A n (t) = P n V{l-P n )e-^ + v)\l-P n )VP n (3.3) 
B n (t) = P n e- l ^+ ir )tp n (34) 

C n (t) = P n V(l - P n )e- i(Ao+ ^ t P„ (3.5) 

and the "residual force" operator F res (t) corresponding to the effect of the 
initial irrelevant part of the state is 

F res (t) = PV(l-P)e- l ^ + ^ t p irr (0)e t{Ao+V)t P 

- Pe-^+^pirr^e^+^il-P^P (3.6) 

For our choice of initial condition, that the system lies within the n-particle 
subspace at t — 0, the residual force term vanishes and our system of equa- 
tions becomes autonomous. 



The lowest order in perturbation expansion The lowest order in the 
perturbation expansion is obtained by replacing the operator e~ lHt in our 
expressions for the kernel for its free counterpart e~ lHot . 

Let us see what happens to the non-local terms when making this ap- 
proximation. The operator C n vanishes since P commutes with e^ H ° T and 
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hence can be made to act on 1 — P. Hence the non-local term in the master 
equation can be written 

- i f dr (-P n V(l - P n )(e- lHoT Ve tH » T )PnP(t - r)e iH ° T P n 

-P n e- iHoT p(t - T)e iHoT e- iHoT Ve iHT (l - P)VP n ) (3.7) 

Now \\e~' lH ° T p(t — T)e lH ° T — p(t)\\ = 0(g), hence by restricting ourselves to 
the second order in the coupling constant we get a master equation that is 
local in time 



.dp 



-£ = [H , p] + [P n VP n , p] + % (L n (t)p + pLl(t)) + F res (t) 



dt 

Here L n (t) stands for the operator: 

L n (t) = P n V(l - P n ) 



f drVi- 



(3.8) 



(3.9) 



where V(r) stands for the time evolution of the operator V according to the 
free Hamiltonian. 

We now set forth to calculate L for the particular cases of n — 1 and 
n — 2. For this we need an expression for the operator W(t) = /q drV(— r). 
The corresponding functional reads 



W(t)(a*, a) = ^ (W abc a a a b a c + 3a* a W a bc a b a c + 3a* a a* b W ab c a c + a* a a* b a* c W abc ) 

(3.10) ' 



with 

W abc ^ 



It 



j(kixi+k 2 x 2 +k3X3) 



l\ (2^,) 1//2 



p -«(^k 1 +^k 2 +^k 3 )t _ i 
x — -(2n) 3 5(k 1 + k 2 + k 3 ) 



"H^ki + ^k 2 + ^k 3 , 
-i(-kixi+k 2 X2+k3X3) 



(3.11) 



X — ■ ■ r(27r) 3 5(k 1 + k 2 + k 3 ) 

-i(-^ki + ^k 2 + ^k 3 ) 



(3.12) 
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W ab [ TT — r -»(-k 1 x 1 -k 2 x 2 +k 3 x 3 ) 

C J M (2UJU-Y/ 2 



x— ■ ■ r(2vr) 3 5(k 1 + k 2 + k 3 ) (3.13) 

-H-^ki + -^k 2 + ^k 3 ) 
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i (s^) 1 / 



^» e i(kixi+k 2 x 2 +k 3 X3) 
e i(aj kl +w k2 +w k3 )t _ ]_ 



— ; r(27r) 3 5(ki + k 2 + k 3 ) (3.14) 

t{uJ kl +LU k2 +UJ k3 ) 



We can easily calculate 



PiVia^a) = |f ((a*a)V abc a a a b a c + 3a* a V a bc a b a c ) 

tW(a*,a) = ^(^(a*a) 2 + V abc a a a b a c 

+3(a*a)a* a V a bc a b a c + 3a* a a* b V ab c a c ) (3.15) 
and from that with repeated application of equation (2.17) we get 

a) = £ {^V abc W ahc + <(V a cd W cd b + V bcd W cda )a b ) (3.16) 



and 



L 2 (t)(a*,a) = 9 - (-(a*a) 2 V abc W abc 



4 V3 

+ (a*a)a* a (2V a cd W cd b + V bcd W cda )a b 
+2 a:«;(2l/>» + V cde W abe )a c a d ) (3.17) 

Renormalization Since we performed a perturbation expansion, it is in- 
avoidable to encounter divergences. We readily see that the term Z = 
V abc W abc as a contraction of two distribution has a 5-function divergence 

Z{t) = / II (9 : , T , T , , (2^(k! + k 2 + k 3 )5(0) (3.18) 
7 . =1 (2u; ki -z(o; kl + ^k 2 + ^k 3 ) 

These terms can be removed from the expression of the operator L(t) by a 
renormalization of the density matrix. To see this it is sufficient to remark 
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that the divergent part of L is proportional to the projection operator P n 
Hence in the master equation these divergencies are contained in the terms 



These terms can be absorbed through a time dependent renormalization of 
the density matrix (essentially wave function renormalization ||) 



which clearly presents positivity, since the term in the exponent is real. 

This is not the only term giving rise to infinities in the master equation, 
but it is more convenient to deal with the others when bringing the master 
equation into a suitable form. 

Having removed these terms we get the following expressions for the op- 
erators 




(Z(t) + Z*(t))p(t) 



(3.19) 




(3.20) 



Li(t) 



—a*F a h a b 
2 a b 




(3.21) 



with 
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e- ikl ( x - x ')(27r) 3 5(k 1 + k 2 + k 3 ) 



(3.22) 




(3.23) 
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(3.24) 
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4 One particle 

Let us now concentrate to the case of the one particle, and in particular 
try to establish the form to which the master equation (3.8) reduces in the 
non-relativistic regime. 

When restricting to the one-particle subspace the reduced density matrix 
is of the form 

p = a* a p\a a (4.1) 

with p a b essentially corresponding the standard form of the density matrix 
p(x, x') in the position representation. Substituting (4.1), (3.) and (3.) in 
the master equation (3.18) we obtain the evolution equation 

i^p\ = h a cP \ - p\h\ + / (F> c b + p\F\) (4.2) 

It is more convenient to consider the density matrix in the momentum rep- 
resentation 

p(k, k') = J «r'e- ikx+lk ' x 'p(x, x') (4.3) 
in which we obtain the transparent form 

i^p(k, k') = K - uv)p(k, k') + ig 2 [F(k, t) + F*(k', t)} p(k, k') (4.4) 
where the function F is given by 

% r dk' 



F(k,t) 



8uJ k J ^k'^k+k' 



e iu > kt (u k cos(^ k / + u k+k /)t + ijfjy + c^k+kp sin(cj k / + cj k+k /)t) - uj k 

X 2 7 \ V? v^-v) 

The function F is actually divergent, but the divergence (logarithmic) is 
essentially contained in its k = argument. Therefore, we can define the 
renormalized function F ren (k, t) = F(k, t) — F(0, t) and absorb the divergent 
part (F(0,t) + F*(0, t) = Z'(t) in a renormalization of the density matrix 

p(t) -> Pren{t) = exp (g 2 J* Z'(t) d?j p(t) (4.6) 
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The non-relativist ic limit To obtain the non-relativistic limit it is suffi- 
cient to obtain the lowest order term in the expansion of F ren (k, t) in powers 
of bfk that is, take the approximation 

*UM) « ^l*=o^' (4-7) 

since the first order derivative of F(k, t) can easily be shown to vanish due 
to symmetry arguments. Writing 

°~ L |k=o(0=7(0^' + iC(*)^' ( 4 -8) 



dttdki 

as we expect again from symmetry requirements, we obtain the master equa- 
tion 

i^ t P=(l + 2img\(t))[Ho)J]+i2mg 2 {H ,p} (4.9) 

Some comments The form of equation (4.9) is clearly unusual , but the 
interpretation is not difficult. Essentially, the effect of the environment is 
to complexify the Hamiltonian, a fact which relates to the non-conservation 
of particles (it is easy to check that this is the term that fails to preserve 
the trace). The term in the anticommutator is a term exhibiting dissipation 
through coupling to energy. 

Also the fact that we have considered a three-field interaction vertex leads 
to the non-appearance of a potential term in the master equation. This will 
be expected to appear when considering at least the n = 3 level. To obtain 
the non-relativistic potential one has to consider theories with two relevant 
field lines at the vertex, as is the cas of QED. 
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